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Abstract 

We derive the thermodynamic entropy of the maen-field spin model in 

the framework of microcanonical ensemble as a function of the energy and mag- 
netization. Using the theory of large deviations and Rugh's microcanonical 
formalism we obtain the entropy and its derivatives and study the thermody- 
namic properties of (j>^ spin model. The interesting point we found is that like 
0^ model the entropy is a concave function of the energy for all values of the 
magnetization, but is nonconcave as a function of the magnetization for some 
values of the energy. This means that the magnetic susceptibility of the model 
can be negative for fixed values of the energy and magnetization in the mi- 
crocanonical formalism. This leads to the inequivalence of the microcanonical 
and canonical ensembles. It is also shown that this mean-field model, displays 
a first-order phase transition driven by the magnetic field. Finally we compare 
the results of the mean-field and (f)'^ spin models. 

1 Introduction 

The microcanonical and canonical ensembles are the two main probability dis- 
tributions with respect to which the equilibrium properties of statistical me- 
chanical models are calculated. Despite the two ensembles model two different 
physical situations, it is generally believed that the ensembles give equivalent 
results in the thermodynamic limit; i.e., in the limit in which the volume of the 
system tends to infinity. The equivalence of the microcanonical and canonical 
ensembles is most usually explained by saying that although the canonical en- 
semble is not a fixed-mean-energy ensemble like the microcanonical ensemble, it 
must 'converge' to a fixed-mean-energy ensemble in the thermodynamic limit, 
and so must become or must realize a microcanonical ensemble in that limit. 
This convergence can be proved to hold for noninteracting systems such as the 



perfect gas and for a variety of weakly interacting systems. For general systems, 
however, neither is this convergence valid nor is the conclusion true concerning 
ensemble equivalence which this convergence is intended to motivate. In fact, 
in the past three and a half decades, numcroiis statistical models have been dis- 
covered having microcanonical equilibrium properties that cannot be accounted 
for within the framework of the canonical ensemble. For systems with short- 
range interactions, the choice of the statistical ensemble is typically of minor 
importance and could be considered a finitesize effect: differences between, say, 
microcanonical and canonical expectation values are known to vanish in the 
thermodynamic limit of large system size, and the various statistical ensembles 
become equivalent . In the presence of long-range interactions this is in general 
not the case, and microcanonical and canonical approaches can lead to different 
thermodynamic properties even in the infinite-system limit . In the astrophysi- 
cal context, nonequivalence of ensembles and the importance of microcanonical 
calculations have long been known for gravitational systems. Although equiv- 
alence of ensembles had been proven only for short-range interactions, it was 
tacitly assumed by most physicists to hold in general. Therefore it came as a 
surprise to many that equivalence does not necessarily hold for long-range sys- 
tems, in particular in the presence of a discontinuous phase transition. Until 
now, the phenomenon of nonequivalent ensembles has been identified and ana- 
lyzed almost exclusively by determining regions of the mean energy where the 
microcanonical entropy hmction is anomalously nonconcave or by determining 
regions of the mean energy where the heat capacity, calculated microcanoni- 
cally, is negative. The existence of such nonconcave anomalies invalidates yet 
another tacit principle of statistical mechanics which states that the one should 
always be able to express the microcanonical entropy, the basic thermodynamic 
function for the microcanonical ensemble, as the Legendre-Fenchel transform of 
the free energy, the basic thermodynamic function for the canonical ensemble. 
Indeed, if the microcanonical entropy is to be expressed as the Legendre-Fenchel 
transform of the canonical free energy, then the former function must necessar- 
ily be concave on its domain of definition. Hence, if the microcanonical entropy 
has nonconcave regions, then expressing it as a Legendre-Fenchel transform is 
impossible. When this occurs, we say that there is thermodynamic nonequiva- 
lence of ensembles. More recently, nonconcave anomalies in the microcanonical 
entropy as well as negative heat capacities have been observed in models of fluid 
turbulence and models of plasmas, in addition to long-range and mean-field spin 
models, including the mean-field XY model and the mean-field Blume-Emery- 
Griffiths (BEG) model. 

One of the motivations for the study of models is the search for soliton. 

Another reason for interest is that $^ models are the simplest systems with 
continuous variables that exhibit a rich phase diagram, with first- and second- 
order phase transitions and a tricritical point. This structure was observed in 
the study of the quantum mechanics of a single site, three-well potential when 
classical, perturbative and mean-field arguments were used and bubble solutions, 
their relation to the phase transitions and the question of their stability, both 
relativistically and non-relativistically. In recent years the model and its 
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application in different physical systems including the following problems have 
been studied extensively: the crossover from a quantum (jfi theory to a renor- 
malized two-dimensional classical nonlinear sigma model, alpha matter on a 
lattice, first-order clcctroweak phase transition (EWPT) due to a dimension-six 
operator in the effective Higgs potential, first-order phase transitions in con- 
fined systems, effective Potential and spontaneous symmetry breaking in the 
noncommutativc <1>^ model, bubble dynamics in quantum phase transitions, the 
canonical transformation and duality in the theory, hermitian matrix model 
for 2D quantum gravity, phase structure of the generalized two dimensional 
Yang-Mills theories on sphere, tricritical Ising model near criticality, sponta- 
neous symmetry breaking at two Loop in 3D massless scalar electrodynamics, 
Ising model in the ferromagnetic phase, statistical mechanics of nonlinear coher- 
ent structures, kinks in the $^ model, growth kinetics in the <I>^ N-Component 
model, stability of Q-balls, the liquid states of pion condensate and hot pion 
matter, instantons and conformal holography, first-order phase transitions in 
superconducting films, field-theoretic description of ionic crystallization in the 
restricted primitive mode. 

This increasing interest in model and its many application in physics, moti- 
vated us to study the statistical mechanics of spin models. 



2 The mean-field $ model and its thermody- 
namic 

The Hamiltonian of the mean- field model is given by the following expression 

1=1 »,j=i 

where qi and Pi are the canonical coordinates of unit mass,moving on a line 
{qi,Pi). The entropy of the system in terms of its mean energy and magnetiza- 
tion is defined as^ 

S(e,m) = lim In / Sieix) — e)6(m(x) — m)dx (2) 

N->-oo N J 

As is sated in [1], if S{€,m) were concave,then one would calculate this 
function from the point of view of the canonical ensemble using the following 

steps 

• Calculate the partition function : 

Z{I3, 7])= (f exp [-l3H{x) - r]M{x)]dx (3) 



where M = Nm{x) 
^In this paper we choose fcg = 1 
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• Calculate the thermodynamic potential defined by 

^(/3,77) = - lim ^\nZ{l3,n) (4) 

where r)) = j3F{f3, 77), where F is the free energy of the system. 

• Obtain S{t, m) by taking the legendre transform of free energy function 
</f(^,??); 

5(e,m) = /3e + r?m-(^(A77) (5) 
with P and ry determined by the equations 

Q Q 

— <f{0,r]) = e , —<f{p,r])=m (6) 

As mentioned earlier these steps are valid only when the entropy is concave. For 
nonconcave entropies we use two methods i.e., a) method of large deviation, b) 
Rugh's method, and compare the results. 

3 Calculation of the entropy using large devia- 
tion method 

We are looking for a set of macro variables or mean-fields such that the energy 
per particle e(x) and the mean magnetization m{x) can be considered as a 
function of these variables: 

e{x) = e{^{x)) , m{x) = fh{^{x)) (7) 

Here we have used iJ,{x) to denote the macro variables collectively. Corre- 
spondingly the entropy can be expressed as a function of these macro variables 
or mean fields: 

Sin) = lim j^^^j ^ifJ'ix) - IJ)dx (8) 

In this method the entropy S{e, m) can be obtained by solving the following 
constrained maximization problem: 

S{e,m)= _ sup S{iJi) (9) 

If S{ii) is a concave function of the given macrostate M, then S{ij) is the 
legendre transform of free energy function defined by : 

Z{\) = j exp [-NX ■ ii{x)]dx (10) 

If S'(/i) is concave, one can find it using the so called macrostate generalization 
of the legendre transform presented in Eqs. (5) and (6) 

s{^i) = ird{\^l-^{X)} (11) 
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Now, let us return to the mean-field (j)^ model. An appropriate choice of 
macrostate is the vector M{m,k^i'), where m, k and v are the mean magneti- 
zation, the mean kinetic energy and the mean potential energy, respectively: 

^ = ^E(^') (12) 



-^E(| + ^) (13) 



AT 6 4 
So we can write as follows: 

'^^^) = k + v-— (14) 

To calculate the entropy we use the legendre transformation. It is shown 

in [2] that the 5(e, m) is given by the following expression: 

fl I ~ 1 

^(e, m) = sup <^ - ln[e -v + —] -\- - ln47re + S{m, u) \ (15) 

4 Rugh's formalism for the mean-field model 

Rugh's microcanonical formalism has been discussed in[l,2]. We assume the 
Hamiltonian H(x;M) in which M is the total magnetization. The microcanonical 
entropy of the system is given by: 

S{E,M) = hi j 5{H{x;M)- E)dx (16) 

where E is the total energy of the system. In Rugh's method the thermodynamic 
quantities are defined through derivatives of the entropy [1,2], and calculated 
by introducing a vector Y in F in such a way that, Y ■ VH = 1, so we have: 

A5(E,M) = (v.y),,^ = ^^ (17) 

and 

^5(E,M) = -/v.(|^^r)\ (18) 



5^ \ 



In general for an arbitrary observable A(x), we have [1]: 

^ {^)e,m = (V.(^l^))s,M - T{E,M) ^^^^'^ ^^^^ 



and 



(20) 
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where (A) denotes the average of the abserable A. By incorporating the mag- 
netization constraint into the Hamiltonian H(x) one can obtain the Hamil- 
tonian H(x;M). This can be done by eUminating one of the qi, for example 
qn = M — 1i- Fo'^ the mean-field model we have : 

N-l 2 1 1,N-1 N-1 2 4 1 

N-l N-l N-l 

i=i i=i i=i 

(22) 

Following Rugh[l], we choose the vector Y as: Y — ^^{pi, ...,pn_i,0, ...,0) 
where Kc is the kinetic part of H(x;M). From Eq.(17) we have: 

1 (23) 

T{E,M) \ 2K I EM 

and from Eq.(18) we have: 

Here ms = J2f=i if "^5 = 7f Si^i if- We define the effective mean-field 
h{E, M) as follows: 

h{E,M) = -T{E,M f^f^^^ (25) 

using Eq.(24), we have: 

/i(£;,M) = -m + r(£,M)/(m3-m5)^^^\ (26) 

magnetic susceptibility and the heat capacity can be calculated by differentiating 
of h{E,M) and T{E,M), see Eqs.(19) and (20). 
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Figure 1: Entropy as a function of the mean magnetization m for different values 
of the mean energy e=0.16, 0.08, 0.04, 0.00 and -0.04 (from top to bottom). 

5 Results of the large deviation method for mean- 
field (f)^ model 

The entropy density ^(e, m) of the mean-field (f)^ model has been shown as a 
function of mean energy e and mean magnetization m for five different values 
of e in fig.l. One of them is above the critical value ec = 0.102 and the rest are 
below. The effective magnetic field h{E, M) has been plotted as a function of 
m for different values of e, in fig. 2. It is worth mentioning that the data in fig.l 
and fig. 2 are obtained using the large deviation method i.e., Eq.(15). In fig. 3 
the magnetic susceptibility has been plotted as a function of magnetization m 
for different values of e. 
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Figure 2: Effective magnetic field as a function of magnetization m for different 
values of e=0.16, 0.08, 0.04, 0.00 and -0.04 (from top to bottom). 

6 Ensemble inequivalence for mean-field 0^ model 

Now, we study the inequivalence of microcanonical and canonical ensemble, for 
the case of (p^ model. This is in close relation with nonconcavity of the entropy 
S{e,m). The effective magnetic field h{E,M) is shown in fig. 2. We observe 
that h{E,M) can be negative for positive values of m when e < Cc. In the 
canonical ensemble, the magnetization m(/3, h) and the magnetic field h have 
the same sign. This means that, in the canonical ensemble, the magnetization of 
the mean-field 4>^ model is always in the direction of the magnetic field. But in 
the microcanonical ensemble the situation is different. h(E,M) and m can have 
opposite signs, so the canonical and microcanonical ensembles are inequivalent. 
For more understanding of the this problem in microcanonical ensemble and its 
relation with nonconcavity of the entropy 5(e, to), we have plotted the magnetic 
susceptibility and to derivative of the entropy in figs. (3) and (4) respectively. 
Since T{E, M) is directly proportional to the kinetic energy, it is always positive 
and this means that the sign of h(E,M) is always opposite to the sign of ^^g^^ ■ 
Then h(E,M) is negative for to > 0. By comparing figs. (3) and (4), we find that 
this happens when the entropy S{e, to) is a concave function of the magnetiza- 
tion e.g., the case for which e = 0.16, then ^"^g^™-* is necessarily negative when 
TO > 0(fig.4), which implies that h{E,M) is necessarily positive (for to > 0). 
This is in agreement with the results of canonical ensemble. Let us now study 
the magnetic susceptibility in the microcanonical and canonical ensemble. In 
the canonical ensemble, the magnetic susceptibility at the constant temperature 
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Figure 3: Magnetic susceptibility as a function of magnetization m calculated 
using the large deviation methode for differennt values of e=0.16, 0.08, 0.04, 
0.00 and -0.04 (from top to bottom). 
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Figure 4: for diffrent values of e=0.16, 0.08, 0.04, 0.00 and -0.04 (from 

top to bottom), culculated using the large divation method. 
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is defined by the following expression: 



5m(/3, h) 
dh 



(27) 



It is easy to show that x^(j3,h) is always positive. In the microcanonical en- 
semble is a function of e and m, and given by: 



dh{e, m) . 



(28) 



where h{e,m) is given by Eq.(25). Using Eq.(25) and (28) one can show that: 



X^(e,m) = -T ^(e,m) 



5e2 



dm? de^ dmde 



(29) 



For the mean-field (f>^ model, is always negative, so if for some values 
of m (when e < ec, sec fig. 3), is positive, then will be negative as can 
be checked by Eq.(29). But we already shown that in the canonical ensmble 

is always positive, this again means that microcanonical and canonical en- 
sembles arc noncquivalcnt. It is worth to mention that this is a consequence of 
nonconcavity of the entropy S{e,m) . 



7 conclusion 

The entropy of the mean-field (p^ model is concave as a function of the energy, 
but is nonconcave as a function of the energy and magnetization. This leads 
to a very important difference between the thermodynamic properties of this 
model in the microcanonical and canonical ensembles. We shown that the ef- 
fective magnetic field in the microcanonical ensemble can have a sign opposite 
to that of the magnetization m, which is in contrast to the case of canonical 
ensemble for which the magnetization m is always in the direction of the applied 
magnetic field. The magnetic susceptibility which in microcanonical ensemble 
is a function of the energy and magnetization can be negative but it is always 
positive in the canonical ensemble. These are two important differences between 
microcanonical and canonical ensemble which make them nonequivalent. The 
mean-field (j>^ model like 0** model displays a first-order phase transition driven 
by the magnetic field in the canonical ensemble which is a consequence of the 
nonconcavity of the entropy, as a function of magnetization (for certain values 
of the energy ) . In addition to display a first-order phase trasition driven by the 
magnetic field the mean-field (j)^ model displays a second order phase transition 
driven by the temperature or the energy. 
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